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Abstract 

We study the B meson spectroscopy allowing the mixture of conventional P wave quark-antiquark 
states and four-quark components. A similar picture was used to describe the new Dj and D s j 
open charm mesons. The four-quark components shift the masses of some positive parity B s j states 
below their corresponding isospin preserving two-meson threshold and therefore they are expected 
to be narrow. Electromagnetic decay widths are analyzed. 
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Heavy-light mesons play in QCD a similar role as the hydrogen atom in QED. This 
analogy provides a simple way to make predictions for their excited states. In the limit Mq — > 
oo heavy-light mesons can be characterized by the spin of the heavy quark, Sq, the total 
angular momentum of the light quark, j q = S q + L, and the total angular momentum, J = 
SQ+jq- For P wave excited states there appear two degenerate doublets: one corresponding 
to j q = 1/2, and the other to j q = 3/2, with quantum numbers J p = + , 1 + and J p = 1 + , 2 + , 
respectively. Those states with j q = 1/2 can only decay via an S wave transition, whereas 
the j q = 3/2 states undergo a D wave transition. Therefore the decay widths are expected 
to be much broader for j q = 1/2 than j q = 3/2 states. We denote the J p = (0 + , l + )j q= i/2 
states as (Aq, AJ and the J p = (1 + , 2 + )j q=3 / 2 states as (Ai, A;), with A = B (for bn states, 
where n stands for a light word quark) or B s (for bs states). Based on this description the 
expected decay properties of the P wave B s j mesons are rather simple, they are summarized 
in Table [H 

A few predictions of masses and widths are available from theory [J @, |, i, H, Q and 
lattice simulations 0], see Table HT1 Similar masses and widths of the L = 1 excited state are 
predicted, around 5.7 GeV for B and 5.8 GeV for B s , except for Refs. ^ and [ij reporting 
lower masses for the B* and B sl . These works use model parameters fitted on the D s meson 
sector and, therefore, they incorporate the anomalies we will comment below. 

From the experimental point of view the spectroscopy of excited mesons containing bot- 
tom quarks is still not well known. Only the ground 0~ and the excited 1~ states are 
well established in the PDG j8[. There are some experimental data by the L3 Collabora- 
tion [9j reporting the first measurement of the B[ and B 2 masses, 5670±10 ± 13 MeV and 
5768±5 ± 6 MeV, respectively. Recently, DO and CDF Collaborations have reported results 
on the spectroscopy of orbitally excited bottom mesons [hJ. CDF found two states, B\ and 
B*, with masses M(B l ) = 5734±3±2 MeV and M(B*) = 5738 ±6 ± 1 MeV. DO also found 
the same states but with slightly different masses, M(Bi) = 5720.8 ± 2.5 ± 5.3 MeV and 
M(B*) - M(B 1 ) = 25.2 ± 3.0 ± 1.1 MeV. In the strange sector CDF reported two narrow 
B sl and B* 2 states with masses M(B sl ) = 5829.4 MeV and M(B* 2 ) = 5839 MeV while DO 
measured only the B* 2 , with a mass of 5839.1 ± 1.4 ± 1.5 MeV. 

The poorly known experimental situation in the open-beauty sector is far from the one 
observed for open-charm mesons. Since 2003, when the _D*j (23 17) and the D s j(2460) were 



discovered by BABAR Collaboration [ljjj, eight cs new states have been reported, more 
than needed to fill the L = doublet and the four L = 1 states. Also the number of cn 
states has grown with the Belle observation [12| of a broad scalar resonance, D$, with a 
mass of 2308 ± 36 MeV and a width T = 276 ± 66 MeV. Some of these new states present 
unexpected masses, quite different from those predicted by quark potential models if a pure 
cq configuration is considered. If they would correspond to standard P wave mesons, their 
masses would be similar to the already known L = 1 cq states, namely around 2.4 GeV 
for cn and 2.5 GeV for cs. They would be therefore above the Dn, D*n and DK, D*K 
thresholds, respectively, being broad resonances. However, the _D*j(2317) and D s j(2460) 
states are below the DK and D*K thresholds and therefore they can only decay through 
the isospin forbidden channels D s tt° and D*it with a very narrow width. In the case of 
the _Dq(2308) the large width observed would be theoretically expected although not its low 
mass. 

These unexpected properties of the _D*j(2317) and the _D s j(2460) mesons have been ex- 
plained in Ref. [l3j assuming a mixture of qq (L — 1, S — 1, 0) and qqqq (L = 0, S = 1, 0) 
components. The reason why four-quark configurations, which mimic quark loop contribu- 



2 



tions, are important in this case is that whereas the qq pair is in a L = 1 state, the four-quark 
state is an S wave and thus its contribution may be relevant. It is well known that the exis- 
tence of four-quark configurations is favored in the heavy-quark sector because, due to the 
coulombic character of the systems, the binding energy augments proportionally to the mass 
whereas the kinetic energy contribution gets reduced when the mass increases. Therefore, 
if sizeable contributions of four-quark structures appear in the D meson sector, they should 
also be present for B mesons. Therefore, one may wonder which are the consequences of the 
mixing between two- and four-quark components in the B meson spectra. 



In the present work we have extended the analysis done on Ref. [13| for open-charm 
mesons to excited P wave open-beauty mesons. Let us first briefly resume the basic features 
of the constituent quark model used [141 ] . The model includes a dynamical quark mass ap- 
pearing as a consequence of the spontaneous breaking of the original QCD SU(3)l<S>SU(3)r 
chiral symmetry at some momentum scale. Once the dynamical quark mass is generated, 
whatever the mechanism, such quarks have inevitably to interact through Goldstone modes 
such that the rotation of the quark fields can be compensated by the boson fields. A simple 



lagrangian invariant under the chiral transformation has been derived in Ref. [15 



L = ip{id- M(q 2 )U^)ip, (1) 

where W 15 = exp(i7r a A a 75// 7r ), 7r a denotes nine pseudoscalar fields (r] Q n, K iy rjs) with 
i =1,...,4 and M(q 2 ) is the constituent mass. The constituent quark mass can be explicitly 
obtained from the theory and parametrized as M(q 2 ) = m q F(q 2 ) with 



F(q 2 



A 2 



(2) 



A 2 + q 2 _ 

Once this has been done U 15 can be expanded in terms of boson fields, 

t/ 75 = 1 + W\ a K a - + (3) 

Jit ^Jjt 

The first term of the expansion generates the constituent quark mass while the second 
gives rise to a one-boson exchange interaction between quarks. The main contribution of 
the third term comes from the two-pion exchange which has been simulated by a scalar 
exchange potential. 

For higher momenta we assume that quarks still interact through gluon exchanges. Fol- 
lowing Ref. [l6| the gluon exchange can be described as an effective interaction between 
constituent quarks given by 

L gqq = i^a^^Xsl) • (4) 

Finally, lattice calculations in the quenched approximation derived, for heavy quarks, a 
confining interaction linearly dependent on the interquark distance. The consideration of sea 
quarks apart from valence quarks (unquenched approximation) suggests a screening effect 
on the potential when increasing the interquark distance [13]. In terms of a potential, the 
color screening produces a gradual decrease of the potential slope, i.e., 

VcoNifij) = {~a c (1 - e-w) + A}(\°i ■ X c j) (5) 

where A is a global constant to fit the origin of energies. At short distances this potential 
presents a linear behavior while it becomes constant at large distances. 
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Using this model we have solved the Schrodinger equation for the two- and four-body 
problems. The two-body case has been solved exactly, while to solve the four-body case 
we have used a variational method with a radial trial wave function taken as the most 
general combination of generalized gaussians [18|. This wave function includes all possible 
flavor-spin-color channels contributing to a given configuration 19| . 



To describe the mixing between qq and four-quark components we use the simplest version 
of the coupled channel model. Namely, it is assumed that a meson state is given by 

W) = 12 a i + 12 Pi ( 6 ) 

i 3 

where q stands for quark degrees of freedom and the coefficients and /3j take into account 
the mixing. The meson systems could then be described in terms of a hamiltonian H = 
Hq + Hx, being 

H °=( H <? H ) andtf 1= ( ° V ^mm\ (7) 

V u n qqqqJ \ ' qq^qqqq u / 

where H is the constituent quark model hamiltonian described above and Hi takes into 
account the mixing between qq and qqqq configurations. It includes the annihilation operator 
of a quark-antiquark pair into the vacuum, that can be obtained from the 3 Pq model. Since 
this model depends on a vertex parameter, we determine this parameter by looking to the 
quark pair that is annihilated (and not to the spectator quarks that will form the final 
qq state). Therefore we have taken V q q^ qq qq = 7. If this coupling is weak enough one 
can solve independently the eigenproblem for the hamiltonians H qq and H qq qq, treating Hi 
perturbatively. To ensure that the perturbative treatment is justified, 7 cannot take all 
possible values, being restricted to \^/{Ej P — Ejp 1 )] 2 < 1. This restriction will limit the 
energy range of the mixed states once the unmixed energies are calculated. The parameter 
7 has been fixed in Ref. [ljjj to reproduce the mass of the _D*j(2317), being 7=240 MeV. 

Let us first compare the experimental results with the predictions of our model for the bn 
and bs quark pairs. In Table II III it can be seen that the model nicely reproduces the known 
experimental data. The two states which have not been measured, (Aq, AJ lie above the 
BK = 5774 MeV and B*K = 5820 MeV thresholds for the strange sector, and above the 
Bit = 5417 MeV threshold for the non-strange one. 

Once the mixing of the qq pairs with the bqqq states is considered, the J p = + and 
1 + states acquire almost a 30% of four-quark component (see Table HV|) . Without being 
dominant, this component is the responsible for shifting the mass of the unmixed states 
below the BK and B*K thresholds. As a consequence, the only allowed strong decay to 
B*ir would not preserve isospin and the resonances are expected to be narrow, of the order 
of a few keV. A second bs J = 1 + resonance appears at M = 5857 MeV, with almost 99% 
of qq component which may correspond with the new B s i state reported by CDF with a 
mass around 5829 MeV. The fourth state appears at 6174 MeV. A similar calculation in 
the non-strange sector is much more involved and we can only predict the existence of a 
Bq state with M = 5615 MeV and 48% of four-quark component and 51% of bn pair. The 
lowest state, representing the Bq, is 200 MeV above the isospin preserving threshold Btt, 
therefore it would be broad. The orthogonal state appears higher in energy, at 6086 MeV, 
also with an important four-quark component. 

Our results are in agreement with Refs. H and [i| obtained with different approaches. 
Ref. 0] interprets the two D s j resonances as the chiral partners of the (0~, 1~) ground state 
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spin multiplet. Fitting the chiral mass gap between the (0 _ , l - ) and (0 + , 1 + ) multiplets to 
the -D*j(2317) experimental mass, they reproduce the -D*j(2460) mass. Assuming that in 
the Mq — > oo limit the chiral mass gap for the charm and bottom sectors is the same, they 
predict the masses for the B meson chiral multiplet (0 + , 1 + ). In Ref. [ij the authors impose 
invariance under heavy quark spin-flavor and chiral transformations to build an effective 
QCD lagrangian. The effective parameters are determined using as inputs the experimental 
D, D s and B s masses together with the assumption that the mass splitting between positive 
and negative parity doublets is the same in the charm and bottom sectors. Their results are 
summarized in Table HT1 These calculations strengthen the idea that the anomalies observed 
in the charm sector must appear in the bottom one. 

Our approach is similar to Ref. j20| considering one loop chiral corrections to calculate 
the _D*j(2317) mass. Such corrections lower the mass of the + state and therefore can ac- 
count for the unusually small mass of the -D*j(2317) and the small mass difference between 
_D*j(2317) and _Dq(2308). However, our approach has the additional property that consid- 
ering four-quark configurations, besides playing a similar role to the one loop corrections, 
it augments the number of states. This provides a plausible interpretation of the recently 
measured -D*(2860) [2l[ as the orthogonal partner of the mixture of conventional P wave 
quark- ant iquark states and four-quark components describing the _D*j (2317) (22| . 

In the L3 Collaboration results 0] the masses of the Bq and B\ mesons, although not 
explicitly given, are constrained by heavy quark effective theory relations [23J, M{B^) — 
M(.Bi) « M(B[) - M(Bl) pa 12 MeV. This allows to estimate the masses of the B* and 
B\ mesons from the B[ and B\ experimental masses, obtaining M(B^) pa 5658 MeV and 
M(B{) pa 5756 MeV. The B\ mass value agrees with the recent measurement of CDF and 
DO Collaborations l(j and with our prediction (see Table IHip whereas the -Bq mass is very 
close to our I = 1/2, J p = + state with a 48% of four-quark component (see Table HVl) . 
Furthermore the L3 Collaboration observed an excess of events above the expected back- 
ground in the 5.9-6.0 GeV region of the Bit mass spectrum, what might correspond to our 
second / = 1/2, J p = 0+ mixed state at 6086 MeV. The L3 data seem to indicate that the 
L — 1 excited B mesons show the same behavior as the corresponding excited states in the 
open charm sector. 

The isovector four-quark states do not couple to the bs system and they are therefore 
much higher in energy, being both close to 6 GeV and above the strong decay threshold. 
Although they should be very broad, and therefore very difficult to observe, any indication 
of its existence would be a definite signal of the presence of four-quark states in the heavy- 
meson sector. 

The structure of the excited B mesons can also be analyzed studying their electromagnetic 
decay widths. The formalism needed to evaluate the T[bq — > bq + 7], T[bnqn — > bnqn + 7], 
and T[bnqn — > bq + 7] widths has been described in Ref. [131 ] . We compare in Table |V] 
our results for the radiative transitions of the + and 1 + states for the cases of a pure 
bq structure, QM(bq), a mixed one, QM(bq + bqqq), and those of Ref. [3J. The difference 
between the QM(6g) results and those of Ref. [3j are due to phase space. The difference 
of the previous two cases with the QM(bq + bqqq) results can be traced back to the more 
involved wave function of this mixed case. The mixing among the 3 P(bq), 1 P(bq) and the 
four-quark components generates a very small 3 P(bq) probability for some particular states, 
as for example the J p = 1 + , of the order of 1%. The electromagnetic decay of the 1 P(bq) 
component to the J p = 1~ state is forbidden. The transition from the four-quark component 
to the meson-photon state does only occur for a very particular component of the tetraquark 
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wavefunction: the one where the light quark- ant iquark pair is in a color singlet, spin one, 
isospin zero state, with a very small probability for the state under consideration. The 
most significant consequence is the suppression predicted for the 1 + — > 1~ + 7 decay as 
compared to the 1 + — > CT + 7. A ratio ~ 1 has been obtained in Ref. 0] and in 

our results with a pure bq structure. For the mixed case we find a much larger value for this 
ratio, j+I^i-^ ~ 100 due to the small l 3 Pi cs probability of the 1 + state. In view of these 
predictions, once experimentally measured, the electromagnetic decay widths would be an 
important diagnostic tool to understand the nature of these states. 

In brief, we have performed an exploratory study of the L = 1 excited B mesons in terms 
of two- and four-quark components based in our experience on the open-charm mesons. Our 
results agree with the recently measured B meson states by CDF and DO Collaborations. 
In addition we predict the existence of two resonances, B* s0 and B' sl , with almost 30% of 
four-quark component, which lie below the BK and B*K thresholds, respectively. Thus, the 
only allowed strong decays would violate isospin and the resonances would be narrow. In the 
non-strange sector we did not find such narrow resonances but our results give support to the 
L3 Collaboration findings. Therefore the mixing between two and four-quark components, 
which explains the unexpected low masses of _D*j(2317) and -D s j(2460) open-charm states, 
would also play a relevant role in the open-beauty sector. We also found that the ratio 
prz^p+^ would provide an experimental signature of the proposed structure. 

We encourage experimentalists on the measurement of the spectroscopic and electromag- 
netic properties of the positive parity Bj and B sJ states, that would help to clarify the 
exciting situation of the open-bottom and open-charm mesons. 
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TABLE I: Expected decay properties of the P wave B s mesons. J is the total angular momentum, 
P its parity, and j q the total angular momentum of the light quark. 



j q J p 


State Decay mode Width 


1/2 0+ 
1/2 1+ 
3/2 1+ 
3/2 2+ 


B* BK broad 
B' sl B*K broad 
B s i B*K narrow 
B* 2 BK, B*K narrow 



TABLE II: bn and bs spectra, in MeV, using different approaches. 



bn 


Rcf. [1] Rcf. [2] Ref. [3] Ref. [4] Ref. [5] Rcf. [6] 


Ref. [7] 


B* Q 


5738 


5760 


5627 


5700 








B[ 


5719 


5780 


5674 


5750 








Bi 


5757 


5780 




5774 


5780 


5755 




B* 2 


5733 


5800 




5790 


5846 


5767 




bs 


Rcf. [1] Ref. [2] Ref. [3] Ref. [4] Ref. [5] Ref. [6] 


Ref. [7] 


B s o 


5841 


5830 


5718 


5710 






5756±31 


B' sl 


5831 


5786 


5765 


5770 






5804±31 


B s i 


5859 


5786 




5877 


5886 


5834 


5892±52 


B s2 


5844 


5808 




5893 


5899 


5846 


5904±52 
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TABLE III: bs and bn quark model (QM) masses, in MeV. Experimental data (Exp.) are from 
Ref. M, except for the states denoted by a dagger from Ref. [lOf] and by a double dagger from 
Ref. ff. 



nL J p 


QM (6s) 


Exp. 


QM (bn) 


Exp. 


is o- 


5355 


5369.6±2.4 


5281 


5279.2±0.5 


1S l- 


5400 


5416.6±3.5 


5321 


5325.0±0.6 


1P 0+ 


5838 




5848 




IP 1+ 


5837 1 
5869 J 


5829.4±0.2 ±0.6+ 


5768 1 
5876 J 


' 5734 ± 3 ± 2+ 

5720.8 ±2.5 ±5.3+ 
k 5670±10±13 tt 


IP 2+ 


5853 


j 5839.6 ± 0.4 ±0,5+ 
[ 5839.1 ± 1.4 ± 1.5+ 


5786 < 


' 5738 ± 6 ± it 

5746.0 ±2.5 ±5.3+ 
k 5768 ± 5 ± 6+t 



TABLE IV: Masses (QM), in MeV, and probability of the different wave function components for 
B s mesons once the mixing between bq and bqqq configurations is considered. 



7 = 


I = 1/2 


J p = 0+ 


J p = 1+ 


J p = 0+ 


QM 


5679 6174 


QM 


5713 5857 


QM 


5615 6086 


P(bnsn) 
P(bs l3p ) 
P(bs 23p ) 


0.30 0.51 
0.69 0.26 
~0.01 0.23 


P(bnsn) 
P(6a x ip) 
P(bs l3p ) 


0.24 ~ 0.01 
0.74 ~ 0.01 
~ 0.01 0.99 


P(bnnn) 
P(bn 1P ) 
P{bn 2P ) 


0.48 0.46 
0.51 0.47 
~ 0.01 0.07 



TABLE V: Comparison of radiative decay widths (keV) assuming only a pure bq structure, QM(bq), 
a combination of two- and four-quark components, QM(bq + bqqq), and those of Ref. 0]. 



Transition 


QM(bq) QM(bq + bqqq) Ref. [3] 


0+ -► 1~ +7 
1+ -> 1- + 7 
1+ -> 0" ±7 


171.4 31.9 58.3 
75.6 0.6 56.9 

106.5 60.7 39.1 


l+^0"+7 
1+-+1-+7 


1.41 101.17 0.69 



s 



